Theory of the Magnetic-Field- Induced Insulator in Neutral Graphene 
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Recent experiments have demonstrated that neutral graphene sheets have an insulating ground state in the 
presence of an external magnetic field. We report on a ;r-band tight-binding-model Hartree-Fock calculation 
which examines the competition between distinct candidate insulating ground states. We conclude that for 
graphene sheets on substrates the ground state is most likely a field-induced spin-density-wave, and that a charge 
density wave state is possible for suspended samples. Neither of these density-wave states support gapless edge 
excitations. 



I. INTRODUCTION 

The magnetic band energy quantization properties of 
graphene sheets lead to quantum Hall effects ''^'^ (QHEs) 



with 



ve^/h at filling factors v = 4(n + l/2) 



■ ■ ■ , —6, —2,2,6, • • • for any integer value of n. The factor of 
4 in this expression accounts for a graphene sheet's two-fold 
valley and spin degeneracies. When Zeeman spin-splitting of 
Landau levels is included, quantum Hall effects are expected 
at the remaining even integer values of V, including the neutral 
graphene v = case. The v = quantum Hall effect of neu- 
tral graphene systems is interesting from two-different points 
of view. First of all, the transport phenomenology of the quan- 
tum Hall effect^ is different at v = because of the possi- 
ble absence of edge states. Indeed the initial experimental 
indications^ that a v = quantum Hall effect occurs in neutral 
graphene did not exhibit either the clear pleateau in p^, or the 
deep minimum in Pvx which are normally characteristic of the 
QHE. Secondly, although a quantum Hall effect is expected 
at V = even for non-interacting electrons, the large energy 
gaps identified experimentally suggest that interactions play 
a substantial role in practice. Gaps due entirely to electron- 
electron interactions in ordered states are in fact common^i^ 
in quantum Hall systems when two or more Landau levels are 
degenerate. Partly for this reason, a number of different sce- 
narios have been proposed^-^' ^-^i'"'''-''^''^''^''^''^ in which the 
gap at V = is associated with different types of broken sym- 
metry within the four quasi-degenerate Landau levels near the 
Fermi level of a neutral graphene sheet. The prevailing view 
has been that the ground state is spin-polarized, with partial 
filling factors Va equal to 1 and — 1 for majority and minor- 
ity spins respectively. This state has an interesting edge state 
structure identical to that of quantum-spin-Hall systems, and 
transport properties in the quantum Hall regime that are con- 
trolled by the properties of current-carrying spin-resolved chi- 
ral edge-states.**''** 

The simplest picture of strong-field physics in nearly- 
neutral graphene sheets is obtained by using the Dirac- 
equation continuum model and neglecting interaction-induced 
mixing between Landau levels with different principal quan- 
tum number n. In this model, electron-electron interactions 
are valley and spin-dependent. When Zeeman interactions 
and disorder are neglected, the broken symmetry ground state 
consists '2 of two-filled n — Landau levels with arbitrary 
spinors in the 4-dimensional spin-valley space. This family 



of states is favored by electron-electron interactions because 
of Fermi statistics which lowers Coulomb interaction energies 
when the orbital content of electrons in the fermion sea is po- 
larized. When Zeeman coupling is included, it uniquely se- 
lects from this family the state in which both n — valley 
orbitals are occupied for majority-spin states and empty for 
minority-spin states. The interacting system ground state is 
then identical to the non-interacting system ground state, al- 
though interactions are expected'^ to dramatically increase the 
energy gap for charged excitations. 

This argument for the character of the ground state appears 
to be compelling, but its conclusions are nevertheless uncer- 
tain. First of all. Landau-level mixing interactions are nor- 
mally stronger than Zeeman interactions, and could play a 
role'^. In addition, although correction s^'"'^'' to the contin- 
uum model for graphene are known to be small at experi- 
mental field strengths, they could still be more important than 
the Zeeman interactions. Suspicions that the character of the 
ground state could be misrepresented by the n = continuum 
model theory have been heightened recently by the work of 
Ong and collaborators, who found a steep increase in the Dirac 
point resistance^' with magnetic field and evidence for a field- 
induced transition to a strongly insulating state at a finite mag- 
netic field strength. Somewhat less dramatic increases in re- 
sistance at the Dirac point have also been reported by other 
researchers . ^3i24,25 

In this article we attempt to shed light on the ground 
state of neutral graphene in a magnetic field by performing 
self-consistent Hartree-Fock calculations for a TT-orbital tight- 
binding model. In the continuum model, Hartree-Fock the- 
ory is known'^ to yield the correct ground state. By us- 
ing a ;r-orbital tight-binding model we can at the same time 
conveniently account for Landau-level mixing effects and 
systematically account for lattice corrections to the Dirac- 
equation continuum model. As we will discuss at length be- 
low, it is essential to perform the mean-field-theory calcula- 
tions with Coulombic electron-electron interactions, and not 
the Hubbard-like interactions commonly used'"'^^ with lattice 
models. One disadvantage of our approach is that our cal- 
culations are feasible only at magnetic fields strengths which 
are stronger than those available experimentally. We therefore 
carefully examine the dependence on magnetic field strength 
and extrapolate to weaker fields. We conclude that under typi- 
cal experimental conditions the most-likely field-induced state 
of neutral graphene on a Si02 substrate is an spin-density- 
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wave state, and that suspended samples might have a charge 
density wave state. Neither of these orderings support edge 
states in the V = gap. We also discuss the magnetic field 
dependence of different contributions to the total energy and 
estimate a critical value of perpendicular and tilted magnetic 
field at which Zeeman splitting will bring about a phase tran- 
sition to a solution with net spin polarization which does sup- 
port edge states. 

Although our calculation captures some realistic features 
of graphene sheets that are neglected in continuum models, 
it is still not a complete all electron many-body theory. In 
particular we neglect the carbon a and a* orbitals whose po- 
larization is expected to screen the Coulomb interactions at 
short distances. Because of our uncertainty as to the strength 
of this screening, our conclusions cannot be definitive. We 
nevertheless hope that our calculations, in combination with 
experiment, will prove useful in identifying the character of 
the field-induced insulating state in neutral graphene. 

Our paper is organized as follows. In Section |ll]we explain 
in detail the model which we study which has two parame- 
ters, a relative dielectric constant which accounts for the 
dielectric environment of the graphene sheet, and on on-site 
interaction U which accounts for short-distance screening ef- 
fects, for example by a-band polarization. Our main results 
on the competition between different ordered states are pre- 
sented in Section |lll] In Section |lV]we turn to a discussion 
of the electronic structure of neutral graphene ribbons, pay- 
ing particular attention to their edge states which play a key 
role in most quantum Hall transport experiments. Finally in 
Section [V]we summarize our main conclusions. 



II. INTERACTING-ELECTRON LATTICE MODEL FOR 
GRAPHENE SHEETS 

A. Non-Interacting Electron ;r-band model 

We first comment briefly on the TT-band tight-binding model 
of graphene in the presence of a magnetic field i^^i^^'^^-^^-^^ 
Each carbon atom on graphene's honeycomb lattice has three 
near-neighbors with ;r-orbital hopping parameter t = —2.6eV. 
Magnetic field effects are captured by a phase factor in the 
hopping amplitudes: t ^ t x e^""^ where = (e/ch) J Adl 
depends on line integral of the vector potential A along a tra- 
jectory finking the two lattice sites. When the dimensionless 
magnetic flux density is =BShex/ 00 = 1 where q is an in- 
teger and (j)Q = cli je is a magnetic flux quantum, it is possible 
to apply Bloch's theorem in a unit cell which is enlarged by 
a factor of q relative to the honeycomb lattice unit cell. (The 
honeycomb lattice unit cell area = ^/3a^ /2 and a = 2.46A 
for a graphene sheet.) Lattice model Landau levels have a 
small width which increases with magnetic field strength and 
reflects magnetic breakdown effects neglected in the contin- 
uum model. 

The ground state energy density differences discussed be- 
low scale approximately as powers of the magnetic length 
is, defined by 2n£gB = 0o- (^b and q are related by £b = 
{Shexq/^nY^^ = 0.371y/qa = 0.913^ A.) In a continuum 



model description the density contributed by a single full Lan- 
dau level is 1 /2k£\ and the energy of the n''' Landau level is 
given by £„ = ±.2hvF \/\n\/ is where vp ~ ^/3at/2h is the 
Fermi velocity of graphene. All energy levels evolve with 
magnetic field except for the n = level, Eq — 0?^ When the 
n''' Landau level is full it contributes E„/ (2n£g) to the energy 
density. From this we immediately see that in the weak-field 
limit important energies tend to scale as ij^^ °^ B^^^. It is easy 
to show, for example, that the magnetic-field dependence of 
the total band energy of a neutral non-interacting graphene 
sheet is given by E (Ib) — aiii„/£g where a^-,-,, = 2.65 eV -A^. 
This non-analytic field-dependence is responsible for the di- 
vergent weak-field diamagnetic response ({dEfot / dB) /B) of 
graphene discussed some time ago by McClure^^. We show 
below that when interactions are included, the energy differ- 
ences between competing field-induced-insulator states also 
tend to vary as £b^. 



B. ;r-band Model Effective Interactions 

It is clear from previous analysis of lattice-corrections to 
continuum model a^i^'i'^i'^" and from lattice-model calcula- 
tions based on extended Hubbard models'" that conclusions 
on the nature of the field-induced insulating ground state are 
very dependent on the effective electron-electron interactions 
used in a TT-band lattice model of graphene. In particular, 
it seems clear that the long-range \/r Coulomb interaction 
tail is essential. We approximate the interaction between k- 
orbitals located at sites separated by a distance li by V [d) = 
1 / {er\/al + cP-) where Oo = a/ (2\/3), the bonding radius of 
the carbon atoms, accounts approximately for interaction re- 
duction due to ;r-charge smearing on each lattice site,'*" and 
Er accounts for screening due to the dielectric environment 
of the graphene sheet. (Here energies are in Hartree (e^/flg) 
units and lengths are in units of the Bohr radius ob-) The on- 
site repulsive interaction parameter, U, is not well known and 
we take it to be a separate parameter We motivate the range 
of values considered for this interaction parameter below. The 
value chosen for can also represent in part screening by a 
orbitals neglected in our approximation, or be understood as 
an ad-hoc correction for overestimates of exchange interac- 
tions in Hartree-Fock theory. Although we study a range of 
values for this interaction parameter model in order to test the 
robustness of our conclusions, we believe that a value of ^ 4 
is normally appropriate for graphene sheets placed on a dielec- 
tric substrate. For practical reasons we truncate the range of 
Coulomb interaction in real space at li = L„,ax = 6.5a. This 
type of truncation is especially helpful when treating systems 
without periodic boundary conditions and allows us to avoid 
problems due to slowly converging sums in real space that can 
otherwise be treated through the Ewald sum method.^' Trun- 
cation of the Coulomb interaction at a reasonably large L„,(„ 
must however be applied with utmost care in order to obtain 
solutions consistent^^-^^ with the limit 

In considering appropriate values for the on-site interaction 
U we can start from the Coulomb interaction energy at the 
carbon radius length scale which is ^ 20eV, while this es- 
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timate can be reduced if one considers a charge distribution 
corresponding to a /j-orbital. In fact an estimate from the first 
ionization energy and electron affinity gives f/ = 9.6eV.27 It IS 
known that the effective on-site interaction strength is greatly 
reduced from this bare value in the solid state environment be- 
cause of screening by polarization of bound orbitals on nearby 
carbon atoms. We consider values of U between 2eV and 6eV, 
bracketing values deemed appropriate by a variety of different 
researchersi22i22j^. A larger value of U increases the inter- 
action energy cost of any charge-density-wave (CDW) state 
which might occur The direct interaction energy is zero when 
all carbon sites stay neutral, but can be positive or negative in 
CDW states. In the CDW states we discuss below electron 
density 5n is transferred between A and B honeycomb sublat- 
tices. In this state the direct interaction energy per site is 



5Eni = 



U+Y,Vid,j)-Y,Vid,j) 

j&A jeB 



(1) 



where dij is the distance between lattice sites / and j, U = 
V{dii) and / is a fixed label belonging to sublattice A. The 
largest terms in Eq.([Tll are the repulsive on-site interactions 
which are proportional in our model to U and attractive ex- 
citonic interaction between electrons on neighboring opposite 
sublattice sites which are inversely proportional to e, . Us- 
ing an Ewald technique to sum over distant sites we find that 
dEoi is positive for e, • C/ > 13.05eV. (The corresponding cri- 
terion for the truncated Coulomb interactions we use in our 
self-consitent-field calculations is • t/ >= 12.23eV; the dif- 
ference between the right-hand-side of these two equations is 
one indicator of the inaccuracy introduced by truncating the 
Coulomb interaction.) When Er-U < 12.23eV the CDW state 
is stable unless band and exchange energies support a uniform 
density state 

Given the band structure model and the interaction model, 
the Hartree-Fock mean-field-theory calculations for bulk 
graphene sheets with periodic boundary conditions and for 
graphene ribbons reported in the following sections are com- 
pletely standard.^' The band quasiparticles are determined by 
diagonalizing a single-particle Hamiltonian which includes 
direct and exchange interaction terms. The direct and ex- 
change potentials are expressed in terms of the occupied 
quasiparticle states and must be determined self-consistently. 
(We do not quote the detailed expressions for these terms 
here.) Since the Hartree-Fock equations can be derived by 
minimizing the total energy for single Slater determinant 
wavefunctions, every solution we find corresponds to an ex- 
tremum of energy. The iteration procedure is stable only if 
the extremum is a minimum so we can be certain that all the 
solutions found below represent local energy minima among 
single Slater determinant wavefunctions with the same sym- 
metry properties. 



III. FIELD-INDUCED INSULATING GROUND STATES 
A. Identification of Candidate States 

At zero-field band energy favors neutral graphene states 
without broken symmetries, and there is no compelling evi- 
dence from experiment that they are induced by interactions. 
In a perpendicular magnetic field, however, the systems is 
particularly susceptible to the formation of broken symme- 
try ground states because of the presence of a half-filled set 
of four-fold spin (neglecting Zeeman) and valley degenerate 
Landau levels with (essentially) perfectly quenched band en- 
ergy. Although the final ground state selection probably rests 
on considerations that it fails to capture, the n = continuum 
model captures the largest part of the interaction energy and 
most of the qualitative physics. The ground state is formed by 
occupying two of the four n = Landau levels, selected at ran- 
dom from the four-dimensional orbital space, and producing 
a gap for charged excitations. 

Three representative broken symmetry states are illustrated 
in Fig. 1 . Because n — Q Landau levels orbitals associated with 
different valleys are completely localized on different honey- 
comb sublattices, a charge density wave (CDW) solution re- 
sults when n = orbitals are occupied for both spins of one 
valley. (When Landau level mixing is neglected valley indices 
and A or B sublattice indices are equivalent.) This state low- 
ers the translational symmetry of the honeycomb lattice in a 
way which removes inversion symmetry. The other extreme 
is a spontaneously spin-polarized uniform density state (ferro 
- F) in which n — Q orbitals are occupied in both valleys but 
only for one spin-component. A third type of broken symme- 
try state, the spin-density-wave (SDW) state, has both broken 
inversion symmetry and broken spin-rotational invariance. In 
the cartoon version of Fig. 1, n ~ electrons occupy states 
with one spin-orientation on one sublattice and the opposite 
spin-orientation on the other sublattice. Possible broken sym- 
metry states, some at other filling factors, had been discussed 
previously by several authors.^''" These three states are all 
contained within the « = continuum model family of ground 
states whose degeneracy is lifted by by lattice nd Landau level 
mixing effects. In the self-consistent mean-field-theory cal- 
culations described in detail below, the three states identified 
above all appear as energy extrema in our collinear-spin study. 



B. Energy Comparisons 

In order to examine the physics behind the competition be- 
tween the candidate ground states we decompose the total en- 
ergy for all three contributions into band, direct interaction, 
and exchange interaction contributions. We have obtained 
self-consistent solutions for all three states over a range of on- 
site interaction U and the dielectric screening values. Be- 
cause of kinetic-energy quenching in the (essentially degen- 
erate) n Landau level, the interaction strengths required 
to drive the system into an ordered state are essentially zero. 
The key question, then, is which state is favored. In Fig. ( 1) 
we illustrate how the energy differences between the three 
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FIG. 1: (Color online) Upper panel: Schematic representation of valley polarized charge density wave (CDW), spin density wave (SDW) and 
ferromagnetic (F) spin polarized broken symmetry solutions that can be obtained in a self-consistent mean field calculation of graphene under 
a perpendicular magnetic field at half filling. Each arrow represents the filling of one n = Landau level of a given spin and valley. Lower 
panel: From left to right iJ^DW _ j7SDW ^ £CDW _ £F ^SDW _ j^j-gj energy differences per electron in eV as a function of the onsite 
repulsion U and obtained from a data mesh of 9 x 10 points, calculated neglecting the Zeeman term and for a magnetic field of 792 Teslas 
corresponding to 1/100 of a flux quantum per honeycomb hexagon. For smaller values of U the CDW solutions are energetically favored 
whereas for larger values of U the SDW solutions are favored in a wide range of e^. The F solutions are never the lowest in energy. The red 
contour lines indicate degeneracy between two different solutions. 
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FIG. 2: (Color online) One dimensional representation of the dispersionless band structure of a graphene sheet under a strong magnetic field 
B = 440r represented in the momentum coordinate k parallel to the narrower direction of the unit cell. The band-gaps follow the i?'/^ scaling 
law expected from the continuum model. The red color is used to represent up spin while blue is used for down spin. Left Panel: Band structure 
for CDW, SDW and F solutions obtained for U = 5eV and e,. = 4. For these interaction parameters the SDW state has the lowest energy and 
the largest gap at the Fermi level. Right Panel: Band structures for U = 5eV and £r = 2. When the on site repulsion is sufficiently weak the 
energetically favored solution corresponds to the CDW state and this solution then has the largest gap. 
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states depend on the model interaction parameters. The re- 
sults in this figure were obtained for ^ = 100 unit cells per flux 
quantum, which corresponds to perpendicular field strength 
B = 792 Tesla. The unit cells in which we can apply peri- 
odic boundary conditions in this case contain 100 x 2 lattice 
sites. The ^-space integrations in the self-consistent Hartree- 
Fock calculations were performed using a 60 A^-point Brillouin 
zone sampling. The self-consistent field equations were iter- 
ated until the total energies were converged to nine significant 
figures. High accuracy is required because the three states 
are very similar in energy since the ordering occurs primar- 
ily in the n = Landau level, involving only 1 % or so of the 
electrons for this value of q. This accuracy was sufficient to 
evaluate energy differences that typically have three signifi- 
cant figures. 

The first point to notice in these plots of energy differences 
is that the two uniform charge density solutions, the F so- 
lution and the SDW solution, behave similarly. The largest 
contrast therefore is between the CDW solution and the SDW 
and F solutions. Focusing first on the CDW/SDW compari- 
son we notice that the SDW state is favored when U is large 
or e, is large. The crossover occurs near Er ■ U ^ 12 eV, 
very close to the line along which 8Eoi changes sign. The 
fact that the CDW/SDW phase boundary occurs very close to 
this line is expected because of kinetic energy quenching in 
a magnetic field. When the non-uniform density CDW state 
is compared with the uniform-density spin-polarized F state 
the phase boundary moves very close to a larger value of this 
product with • U ranging from ^ 14 to ~ 18eV along the 
phase boundary. Evidently the competition between CDW 
and SDW states is based very closely on the direct interaction 
energy, with additional weaker elements of the competition 
entering when the F state is considered. 

Direct comparison between the uniform density SDW and 
F solutions indicates that the latter is favored only at values 
of U and e, which are outside the range of most likely values. 
As discussed in more detail below, we find that the direct in- 
teraction energy in these two states is identical, and that the 
more negative exchange energy of the SDW state overcomes 
a larger band energy. In this case the main difference between 
the energies of the two states arises from Landau level mixing 
effects. As we explain later. Landau level mixing leads to a 
local spin-polarization which is larger in the SDW state than 
in the F state. 



In Fig. ( 1) we have introduced the main trends in the ener- 
getic competition between CDW, SDW, and F states. How- 
ever, as we have explained, these calculations were under- 
taken at field strengths that exceed those available experimen- 
tally. In the following subsection we demonstrate that the field 
dependence of the energy comparisons is extremely system- 
atic so that extrapolations down to physical field strengths are 
reliable. So far we have also ignored Zeeman coupling which 
favors F states. This coupling can be important and is also 
addressed in the following subsection. 
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FIG. 3: (Color online) Total energy per site differences AE,oi, sep- 
arated into kinetic energy AEi,>,, electrostatic energy AEoj and ex- 
change energy AExi contributions as a function of magnetic length 
£b in lattice constant a = 2.46A units. The total energy differences 
were fitted to a €2/2^^^'^ + €28-^ curve. The fitting parameters are 
listed in Table I. Left panel: Energy differences between F and 
SDW solutions KE^'I^DW ^ _ ^SDW j^^^^ ^g^yjj^ ^^^^ ^j,. 

tained with interaction parameter values U = 5 eV and = 4 for 
which SDW is the lowest energy configuration. The more negative 
values of exchange energy in the SDW state compensates the kinetic 
energy penalty related to the inhomogeneous accumulation of the 
electron wave functions at alternating lattice sites. The electrostatic 
energy differences are zero thanks to the uniform electron density 
for both solutions. Right panel: Same as the previous figure but for 
^F/CDW _ _ £CDW -pjjg interaction parameters in this case are 
U = 5 eV and £r = 2 for which CDW is the lowest energy configu- 
ration. When the onsite repulsion U is small enough that the elec- 
trostatic energy penalty for the inhomogeneous charge distribution 
is small, exchange is the main contribution driving the CDW insta- 
bility. However, in the case illustrated here U is so small that the 
electrostatic part of the hamiltonian does play an important role in 
favoring the CDW state. The energy contributions follow a magnetic 
field decay law that deviates more from B^/^ than in the previous 
case because the on-site interaction U plays an essential role. 



200 



100 





u 


= 5eV, Er = 4 






= 5eV, £,. = 2 






>w 



2000 



1000 



7r/6 



7r/3 



tt/2 



FIG. 4: Tilt angle 9 dependence of the critical magnetic field re- 
quired to induce a transition to the F state. The black solid curve and 
dashed blue curve represent the critical magnetic fields starting from 
the SDW and CDW states respectively. In the CDW curve we ob- 
serve a larger deviation from a simple cos(0) law due to a stronger 
influence of lattice scale physics described by the C2 coefficient in 
equation [3] 
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TABLE I: In the upper table we list the values of q™"'''^"'^' in units of lO^'^eV/r^/^ and c^™'*'/'^^'*') in units of IQ-^^eV /T^, obtained 
from fits to the their energy difference with respect to F states as given in equation l|2} for two different values of the interaction parameters U 
and £,-. The critical field estimates are dependent on the reliability of these fits. Note that a crossover between SDW and CDW states can be 
driven by changes in the dielectric screening environment captured by e,-. The lower table lists Be values at which the F states become critical 
according to Eq.l[3]l. 



C. Field Strength and Zeeman Coupling Dependence 



We now turn our attention to the magnetic field dependence 
of the solutions. For this purpose we found self-consistent 
solutions over a range of magnetic fields for two sets of inter- 
action parameters, U = 5eV and = 4 for which the SDW 
solution has the lowest energy, and U = 5eV and Er = 2 for 
which the CDW solution has the lowest energy. The band 
structures of the different possible solutions for these set of 
parameters are shown in Fig. (2) and the field dependences of 
the three energy differences are plotted in Fig. (3). We see 
that every contribution accurately follows a Z?^/^ °^ Ig^ law 
with small deviations that can be accounted for by allowing 
a term proportional to B^. This is the same field-dependence 
law that we discussed earlier for the case of a non-interacting 
electron system. In the continuum model it is guaranteed in 
neutral graphene when electrons interact via the Coulomb in- 
teractions by the fact that both kinetic and interaction energy 
densities then scale as {length)^^; the magnetic field sim- 
ply provides a scale for measuring density. The fact that we 
find this field dependence simply shows that the condensation 
energies of all three ordered states are driven by continuum 
model physics. This is in agreement with the intuitive pic- 
ture of the interaction energy as the product of the number of 
electrons occupying a Landau level which is directly propor- 
tional to B, multiplied by the Coulomb interaction scale for 
electrons in the n = Landau level which is proportional to 
B^/^. The fact that the differences in energy between the three 
states follows this rule suggests that the most important source 
of differences in energy between these states is Landau-level 
mixing, which should not violate the B^/^ law. Small devia- 
tions from this law are expected because of lattice effects. The 
deviations are stronger in CDW solutions than in the SDW 



solutions because of the charge density inhomogeneity at the 
lattice scale present in the former. 

We can draw two important additional conclusions from the 
B^/^ behavior. First of all, lattice effects are not dominant 
effect at the field strengths for which we are able to perform 
calculations, and should be less important at the weaker fields 
for which experiments are performed because the magnetic 
length Ib will then be even longer compared to the honeycomb 
lattice constant. The difference in energy between the three 
states should mainly vary as B^/^ all the way down to zero 
field, provided only that disorder is negligible. (We discuss 
the role of disorder again in Sec. [V]) Our calculations should 
therefore reliably predict the energetic ordering of the states in 
the experimental field range. The second conclusion we can 
make concerns the importance of Zeeman coupling which we 
have ignored to this point. First of all, Zeeman coupling will 
have a negligible effect on the energies of the SDW and CDW 
states since they have a vanishing spin magnetic susceptibility. 
The energy difference per site between the F state and the two 
density-wave states can be written in the form 

AE = £'(™W'/CDW)_£.F 

= B3/2cg7/^^"''cos(0)3/2 

+ b2. (c(^^"'/^^"')cos(0)2-Czcos(0)) (2) 

where B is the total magnetic field strength and 9 is the field 
tilt angle relative to the graphene plane normal. Factors of 
Bcos(0) in this expression therefore account for the perpen- 
dicular field dependence. The second term in Eq. |2] contain 
the contributions that scale with B^. The factor of Bcos(0) 
which appears in the Zeeman term is present because the spin- 
polarization of the F state is proportional to the Landau level 
degeneracy. The coefficients C3/2 and C2 can be obtained 
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by fitting energy differences obtained from numerical solu- 
tions of the self-consistent field equations, like those plot- 
ted in Fig. 3, and depend on the interaction model parame- 
ters as shown in Table I. The Zeeman coefficient in Eq. |2] is 
Cz = 7.3 • 10^'" eV /T^ is independent of interaction param- 
eters. From the above equation we find that the F state has 
lower energy than the spin-unpolarized states for 



B>Bc{e) 



'-3/2 



cos(0) 



(Q-C2COS(0))' 



(3) 



The fields required to achieve an energetic preference for the 
spin-polarized state are smaller at larger tilt angles because 
the orbital energy has a stronger 9 dependence. 

In table I we show the values of C3 /2 and C2 for SDW and 
CDW configurations favored with respect to F for a set of pa- 
rameters of U and e^. We notice that the coefficients dictating 
the critical field transition to F solutions can be made rela- 
tively small if the parameters are near the crossover boundary 
to F states. It is possible that a SDW or CDW to F transition 
could be induced by varying magnetic field. If a transition 
was observed, most likely by a change in transport proper- 
ties as discussed in the next section, it could provide valuable 
input on the effective interaction parameters of the ;r-orbital 
tight-binding model. 



IV. QUANTUM HALL EDGE STATES IN GRAPHENE 
RIBBONS 



a = 2.46A 




Zigzag ribbon unit cell 



FIG. 5: Representation of the unit cell choices for armchair and 
zigzag edge terminated graphene nanoribbons. 



The quantum Hall effect occurs when a two-dimensional 
electron system has a chemical potential discontinuity (a gap 
for charged excitations) at a density which depends on mag- 
netic field. A gap at a field-dependent density necessarily^"'^' 
implies the presence of chiral edge states that support an equi- 
librium circulating current. The current varies with chemi- 
cal potential at a rate defined by the field-dependence of the 



bulk gap density. Most quantum Hall measurements simply 
reflect the property^' that separate local equilibria are estab- 
lished at opposite edges of a ribbon in systems with a bulk 
energy gap or mobility gap. It is immediately clear therefore 
that the v = quantum Hall effect is special since it is due to 
an energy gap at the neutrality point, i.e. at a density which 
does not depend on magnetic field. The issue of whether or 
not the V = gap and associated phenomena should be re- 
ferred to as an instance of the quantum Hall effect is perhaps 
a delicate one. The v = gap is intimately related to Landau 
quantization and in this sense is comfortably grouped with 
quantum Hall phenomena. This view supports the language 
we use in referring to the v = quantum Hall effect. On the 
other hand, since it occurs at a field-independent density, its 
transport phenomena are more naturally viewed as those of an 
ordinary insulator'^ which just happens to be induced by an 
external magnetic field. 

An exception occurs for the F state which does have edge 
states^—, and can be viewed as having v = 1 for majority 
spins and V = — 1 for minority spins. In the simplest case, 
it has two branches of edge state with opposite chirality for 
opposite spin, much like those of quantum-spin-Hall'^ sys- 
tems. In a Hall bar geometry most transport measurements 
are very strongly sensitive to the presence or absence of edge 
states. In order to address edge state physics directly at v = 
we have extended our study from the bulk graphene to the 
graphene nanoribbon case. Tight-binding model solutions for 
a ribbon in the presence of a magnetic field can be obtained in 
essentially the same way as for bulk graphene, with the sim- 
plification that any magnetic field strength preserves the one- 
dimensional ribbon wavevector A; as a good quantum num- 
ber when the gauge is chosen appropriately. This graphene 
ribbon problem in a magnetic field was studied time ago by 
Wakabayashi et. al.^'' and revisited recently within both tight- 
binding^^'^^ and continuu m^'^^'^^ models in order to provide 
a microscopic assessment of the relationship between Lan- 
dau levels and edge states. The general feature of the ribbon 
band structure in the presence of a magnetic field is that those 
states localized near the edges have dispersive bands, whereas 
those in the flat band region are located mostly in the bulk. In 
the case of zigzag edge termination, edge localized states are 
present even in the absence of a magnetic field^^. In the quan- 
tum Hall regime these states are in the non-dispersive band 
region like the bulk localized states and they do not contribute 
to edge currents, although they can interact with other edge 
localized states. 

Fig. 6 explicitly illustrates how the character of the bulk 
broken symmetry is manifested in ribbon edge state prop- 
erties. Because of practical limitations our calculations are 
restricted to moderately narrow ribbons with widths of or- 
der lOnm. In order to properly reproduce bulk Landau level 
quantization in these narrow systems we have to choose mag- 
netic field strengths strong enough to yield magnetic lengths 
Ib ^ 25nm/(Z?[Tesla])'''^ substantially smaller than the ribbon 
width, i.e. fields stronger than typical experimental fields. On 
the other hand if the magnetic fields are too strong, say £b < a 
the levels will be strongly affected by the lattice and the prop- 
erties of the solutions will substantially depart from the behav- 
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FIG. 6: (Color online) Band structure and spin density distributions in armchair ribbons calculated for C =4 eV and £,■ = 4 for a magnetic 
field corresponding to = 4.24a ( B = 606 Tesla) using a 100 x 2 lattices in the unit cell and 80 k-point sampling in the Brillouin zone. The 
red color is used to represent up spin while blue is used for down spin. Upper row: Armchair ribbon band structures for the F, SDW and 
CDW solutions. Only the F solution is metallic, while the other two configurations are insulating. For this choice of parameters SDW is the 
energetically favored state, but the energies of the other states are similar and their band gaps also have a similar size. Lower row: Spin-^ and 
spin-| density distributions across one of the zigzag rows in the unit cell of an armchair ribbon. The influence of Landau-level mixing on the 
three solutions is apparent in the differences between the three sets of density distributions. In the case of AF solutions Landau-level mixing 
enhances the local spin-density, while in the case of the CDW solution the magnitude of the charge density oscillation is reduced by Landau 
level mixing. The dashed horizontal lines represent the occupation in absence of Landau level mixing. 



ior we should expect at weaker magnetic fields, for which the 
continuum model description is approximately correct. For 
field strengths in the appropriate range, we find the same three 
types of self-consistent field solutions as in the bulk calcula- 
tions, namely F, CDW, and SDW solutions. The band struc- 
tures and spin resolved densities presented in Fig. 6 confirm 
that only the F solution has states in the broken-symmetry in- 
duced gap. Hall-bar transport properties for the F configura- 
tion have been discussed by Abanin et al.^ and Fertig et al— 
from a theoretical point of view. The CDW and SDW state 
electronic structure is insulating, both at the edge and in the 
bulk. The band structures of these two-states are similar even 
though their spin-density profiles are quite distinct. 

The plots of spin-t and spin-| partial densities across the 
ribbons hint at some of the physics which selects between the 
three candidate ordered states. In the truncated n~Q Landau- 
level continuum theory, the F state has one excess occupied 
Landau level for majority spins and one deficiency in occu- 
pied Landau levels for minority spins. We see in Fig. 6, that 
the size of these polarizations is not strongly influenced by the 
Landau level mixing effects included in our lattice calculation. 
The n — Q continuum theory SDW state has the same spin ex- 
cesses and deficiencies, but they have opposite sign on oppo- 
site sublattices. We see in Fig. 6, that these order parameters 
are actually enhanced by Landau level mixing effects; inter- 



Landau level exchange effects polarize lower-energy occupied 
Landau level states so that they enhance the SDW pattern. For 
the CDW state on the other hand, the n—Q Landau level ex- 
cess density on one sublattice is suppressed by Landau-level 
mixing. In this case the direct electrostatic interaction is non- 
zero so that the occupied Landau levels away from the Fermi 
energy are polarized between sublattices in the opposite sense 
of the n = levels. The fact that the SDW state is enhanced 
by Landau level mixing explains why it is favored over the F 
state for the interaction parameters used to construct Fig. 6. 

Finally, we comment on the microscopic electronic struc- 
ture at the edges of an F state. As illustrated in Fig. 6, 
the mean-field electronic structure at the edge contains a do- 
main wall.^^ Because textures in this domain wall can'** carry 
charge, the F state edge is however not necessarily insulating 
in the absence of disorder. It was recently argued that impuri- 
ties with magnetic moments near the sample edges can intro- 
duce spin-flip backscattering potentialsi^ whose effectiveness 
is enhanced by the presence of a domain wall. The high resis- 
tances seen experimentally at v = therefore do not necessar- 
ily prove that the ground state is not an F state. The analysis 
in RefTs* was caiTied out within a Luttinger liquid formalism 
whose parameters depend on the domain wall shape. As il- 
lustrated in Fig. 7, our microscopic domain walls exhibit an 
interesting anisotropy in which inner and outer segments of 
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FIG. 7: (Color online) Quantum hall edge domain-wall wave func- 
tions and energy bands corresponding to an F solution calculated 
for a zigzag graphene nanoribbon with 20 carbon atom pairs in the 
unit cell (W = 42. 6A) under a magnetic field of 6 = 5660 Teslas. 
Left panel: In the upper figure we represent the edge state wave- 
functions for the occupied and unoccupied t-spin and |-spin states 
when hybridization is not included. The lower figure represents 
the coefficients of the hybridized domain wall state. All wavefunc- 
tions are plotted at the Bloch wavevector indicated by arrows on 
the right panel. Right panel: Quantum hall edge state bands with 
collinear spin states represented by black and red symbols and the 
non-collinear spin edge states formed by optimizing the hybridiza- 
tion of up and down spin orbitals represented by the blue symbols. 
We observe a clear single-particle energy gap related to the energy 
gained by forming the non-collinear domain wall structure. 



the wall differ. 



V. SUMMARY AND DISCUSSIONS 

We have carried out a mean-field study of graphene's v = 
ground state which aims to shed light on the character of the 
interaction-induced gap that appears experimentally. The fact 
that the ground state has a charge gap at this filling factor can 
be established essentially rigorously'^ within the continuum 
model often used to describe graphene. When interaction in- 
duced mixing between n = and |n| 7^ Landau levels is ne- 
glected in the continuum model, the family of broken sym- 
metry states is related by arbitrary spin and valley pseudospin 
rotations and includes both fully spin-polarized F and spin 
(SDW) and charge (CDW) wave states. Our effort to deter- 
mine the ground state character when Landau-level mixing 
and lattice effects are included, is motivated by the observa- 
tion of magnetic-field-induced insulating transport properties 
and by the expectation that transport properties in the quantum 
Hall regime should be very different for F, SDW, and CDW 
states. 

The lattice model we study has two phenomenological pa- 
rameters, a relative dielectric constant and an on-site in- 
teraction parameter U . The most appropriate values of both 
parameters are somewhat uncertain. Our two main findings 



are that i) Landau-level mixing effects favor the density-wave 
states over the feiTomagnetic state and ii) The competition 
between CDW and SDW states is sensitive to the relative 
strength of on-site and inter-site electron-electron interactions 
and hence on the product £,. • U . Large values of this prod- 
uct increase the direct mean-field energy cost of CDW order 
and favor the SDW state. Exchange energies are stronger for 
density wave states than for ferromagnetic states because or- 
der within the n — Q level, induces order in the full negative n 
Landau levels in the former case, but not in the latter case. For 
graphene on Si02 and other typical substrates, sensible values 
for Er and U suggest that the field-induced state at n = is a 
SDW state. CDW states could occur in suspended graphene 
samples. 

The atomic value of the on-site interaction term in carbon 
is f/ ~ lOeV, so graphene values should be smaller Most 
of the illustrative calculations we have described use either 
U = AeV or U = 5eV . The dielectric constant is ~ 1 for 
free standing graphene. For a Si02 substrate, with dielectric 
constant ~ 4.5, the effective 2D dielectric constant at a sub- 
strate/vacuum interface is ^ 2.5. Because the Hartree-Fock 
approximation overestimates the strength of exchange inter- 
actions, it can be argued that somewhat larger values of 
are appropriate - perhaps ^ 2 for free standing graphene and 

~ 4 or 5 for graphene on a substrate. Because of these un- 
certainties we view U and as effective parameters whose 
values are somewhat uncertain and have made energy com- 
parisons over a wide range of values. 

We are able to complete our calculations only at magnetic 
field strengths much stronger than those available experimen- 
tally. Partly to verify that the field-dependence is system- 
atic, and partly in an effort to estimate the magnetic field 
strengths necessary for Zeeman energies to drive the system 
into a F state, we have fit the energy differences between F 
and density- wave states to the form AE ~ Ct^ijB^I'^ + C^B^. 
Only the first term can be present for a continuum model 
with Coulomb interactions, and we find that this term is in- 
deed dominant. We find that Zeeman coupling at typical fields 
can change the nature of the state only in the parameter range 
where the crossover from CDW to SDW states occurs. Al- 
though we assume collinear states in analyzing the SDW/F 
competition, we presume that the SDW to F crossover is ac- 
tually a continuous one in which the antiparallel spins on op- 
posite sublattices are smoothly rotated until they are aligned. 
Even if the bulk gap remains open the rotation of the spins will 
bring about a progressive closing of edge state band gaps, until 
the gap is completely closed in the F configuration. The mod- 
ulation of the edge charge gap due to Pauli paramagnetism 
might be detected in edge transport experiments. Since a rel- 
atively strong magnetic field is required to induce the insulat- 
ing state in typical samples, we do not expect that it will be 
easy to produce parallel fields large enough to turn the system 
metallic while maintaining this minimum perpendicular field. 
Nevertheless, a study of the parallel-field dependence of trans- 
port properties is likely to hint at the nature of the underlying 
broken symmetry state of the system. 

All of these results ignore the influence of disorder Experi- 
ments appear to show that the transition to the insulating state 
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occurs at magnetic fields above some critical value that be- 
comes smaller when the sample is cleaner This is expected, 
since disorder favors states without broken symmetries. The 
relationship between the minimum field and the mobility was 
carefully examined some time ago/ but can be crudely de- 
scribed using the following simple argument. Assuming un- 
correlated scatterers and using the Fermi golden rule the mo- 
bility is /X o<: 1 /Vjij where ViUs is the typical energy scale of 
disorder. The crossover occurs when disorder strength equals 
the interaction energy scale V^is = Ui„t °= Br'l'^, therefore the 
critical magnetic field between samples with different mobil- 

1 /3 

ity are related through B'^/Bc — {n/ n') ■ One physical pic- 
ture of the strong disorder limit asserts that current flows along 
domain walls which separate disorder-induced electron-hole 
puddles domain walls along which current carrying states 
with V 7^ can dominate bulk transport and suppress the di- 
vergent resistivity. In graphene on substrates the mobility val- 
ues range between 2000 and 25,000 crn^V ^ ' 5^ ' , and values as 
high as 230,000 cm^V^^s^^ have been achieved in annealed 
suspended samples^. Typical crtical fields in samples on sub- 
strates are in the 20 to 30 Tesla range. From the above argu- 



ment we can expect that the critical magnetic fields in sus- 
pended samples should be roughly 2 to 5 times lower. 

An important goal of our work was to shed light on the 
character of the v = edge states. We examined the elec- 
tronic structure of armchair ribbons with CDW, SDW, and F 
states, finding that edge states in the gap are absent for both 
CDW and SDW solutions. Given that the field-induced insu- 
lating state appears to have an extremely large resistance once 
established, it appears likely to us that the experimental state 
does not have edge states and that it therefore must be a den- 
sity wave, as suggested by these calculations. If so, a study of 
the influence of magnetic field tilting might be able to distin- 
quish between SDW and CDW states. 
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